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O i Abstract 

. A general recipe proposed elsewhere to define, via Noether theorem, the variation of energy for 

'sj" ' a natural field theory is applied to Einstein-Maxwell theory. The electromagnetic field is analysed 

. in the geometric framework of natural bundles. Einstein-Maxwell theory turns then out to be 

natural rather than gauge-natural. As a consequence of this assumption a correction term a la 
^ . Regge-Teitelboim is needed to define the variation of energy, also for the pure electromagnetic 

' part of the Einstein-Maxwell Lagrangian. Integrability conditions for the variational equation 

. which defines the variation of energy are analysed in relation with boundary conditions on 

physical data. As an application the first law of thermodynamics for rigidly rotating horizons is 
obtained. 



Introduction 



In a previous paper |28| a recipe to define, via Noether theorem, the variation of Noether conserved 
quantities in natural field theories was proposed. The main result of the theory there developed 
was to describe, for vacuum General Relativity, the quasilocal stress-energy content of a spatially 
bounded region of spacetime with non-orthogonal boundaries. 
^ ' The purpose of the present paper is to generalise the recipe of [28| to Einstein-Maxwell theory 



■ in order to calculate the Noether quasilocal conserved quantities, still remaining in the general 

case of spacetime regions with non-orthogonal boundaries. As a remarkable application we define 
a first principle of thermodynamics for spacetimes admitting an isolated horizon (in the sense of Q ) . 

Using a geometric approach to field theories in the framework of natural theories, it is possible to 
associate to each infinitesimal symmetry ^ on the base manifold M (spacetime of dimension m) 
a conserved quantity Q{S,). This Noether charge is obtained by integrating the Noether current 
on a compact (m — 1) submanifold S of M with boundary dT, and it can be decomposed into a 
bulk term QbuikiO and a surface term Qsurfacei^,)- The bulk term is obtained by integrating on 
S the reduced Noether current, which is a (m — l)-form vanishing on shell (i.e. along solutions). 
The surface term is the integral over dT, of the superpotential related to the symmetry generated 
by ^ which in turn, is a (m — 2)-form algorithmically calculated starting from the Lagrangian of 
the theory |23]. To obtain physically expected values, in analogy with the original prescription 



suggested by Regge-Teitelboim in |43l (when dealing with the analysis of Hamiltonian boundary 
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terms), we have to correct the Noether conserved quantities by suitably adding boundary terms to 
Q(^); this approach was extensively developed for natural and gauge-natural theories in a global 
and geometric framework, called covariant ADM formalism. [23|. 

Following the same idea which leads to the covariant ADM formalism, we can proceed one step 
further by defining the variation of the corrected conserved quantities through the addition of a 
further suitable boundary term r (which is still present in the covariant ADM definition of the vari- 
ation of the conserved quantity), such that SQ{^) turns out to be a pure (off-shell) bulk term. As a 
particular case, if we consider a Cauchy surface E in spacetime and a vector field ^ transverse to S 
the variation of the Hamiltonian can be defined as the variation of the Noether charge associated 
to ^. In this situation the vector ^ is identified with the generator of time fiow, meaning that the 
parameter "t" generating the fiow of ^ itself is identified with time. This technique allows us to 
obtain the Hamiltonian equations of motion and to define the symplectic structure of the theory 
together with its phase space. Variation of energy is simply defined as the on-shell value of the 
variation of the Hamiltonian. A problem which arises is then to formally integrate, if possible, 
the variation of the conserved quantities to obtain the explicit expression for the Hamiltonian and 
for the energy. This problem is related to the choice of suitable boundary conditions as well as 
to the choice of a background solution. The chosen background can be interpreted as a reference 
solution in the one-parameter family of solutions satisfying the same boundary conditions (namely 
the one-parameter family along which variations are performed). 

We remark that the final formula for the variation of the Hamiltonian is independent on divergence 
terms which might be added to the Lagrangian. This is a mathematical and physical advantage 
since we do no longer have to take care about surface terms in the action functional. The expression 
of the variation of the Hamiltonian is unique for each element of the equivalence class [L] formed by 
all the Lagrangians differing one from the other only for divergence terms. Moreover, the variation 
of the Hamiltonian reproduces the correct equations of motion in the whole phase space of the 
theory and not only in a phase space restricted by boundary conditions such as, for example, a 
suitable fall-off requirement on solutions at spatial infinity. 

In this framework, in fact, the symplectic structure of the field theory is uniquely defined by the 
bulk term, generated by the so called reduced symplectic current. There is no contribution to the 
symplectic structure from the boundary terms, which are instead pushed directly into the very 
definition of the variation of the Hamiltonian. 

The variation of energy, as previously remarked, is defined as the on-shell variation of the Hamil- 
tonian. In this framework different definitions of energy for the same system do in fact arise in 
correspondence to different definitions of the variations of the Hamiltonian (and consequently of 
the variations of energy) . Different variations are related to the different vertical vector fields which 
generate different one-parameter families of solutions starting from the same solution. Different 
one-parameter families correspond to different boundary conditions the solutions have to satisfy. 
The variational equation which defines the variation of energy can be (or can be not) integrated, 
depending on the boundary conditions imposed, when integrable different kinds of energy are de- 
fined. 

This viewpoint is in full accordance with the classical treatment of thermodynamical systems. When 
we impose different boundary conditions on the same thermodynamical system, in fact we perform 
each time a different choice between the intensive or extensive variables for the system (in a sym- 
plectic framework we can say that we are choosing the "control parameters") and correspondingly 
we expect to find different energies for the system. 
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In vacuum General Relativity, as shown in |28|, the definition of quasilocal energy via Noether 
theorem leads to the same results previously obtained by using the Trace-K action functional |1C], 
1 11 1 or by using a Hamiltonian analysis as in Q, |37]. The Hilbert Lagrangian is sensitive to the 
formalism presented here. In presence of non-orthogonal boundaries, in fact, the generalized Regge- 
Teitelboim correction term r (we add to the definition of 6H) is fundamental for the definition 
of the variation of energy of the gravitational field. In this way a definition of quasilocal internal 
energy for spatially bounded gravitating systems can be obtained where the reference background 
is properly taken into account. This is obtained by imposing Dirichlet boundary conditions, i.e. the 
metric at the boundary of the world tube is kept fixed. 

In this paper we apply our general recipe PI to the case of Einstein-Maxwell theory. Electro- 
magnetism is here treated as a natural theory. In natural theories each vector field tangent to the 
spacetime manifold can be naturally lifted to a vector field tangent to the configuration bundle 
and moreover this lifted vector field is a symmetry for the Lagrangian |^] . We remark that in the 
natural lift of spacetime vector fields no gauge freedom remains undetermined (quite different is the 
case in gauge theories!). In this well defined geometric framework the Lie derivative of the fields 
(i.e. the electromagnetic potential) is uniquely defined, has a correct mathematical interpretation 
and it is thence possible to define Noether conserved quantities via Noether theorem. 
Maxwell theory, opposite to other gauge theories, can be treated as a natural theory by using a 
suitable representation of the gauge group U{1). The configuration bundle turns out to be a natural 
U{1) bundle; this in turn implies that the configuration bundle is a trivial U{1) bundle, meaning 
that no magnetic charge is allowed: see pl[| . 

In this particular framework the electromagnetic potential is a geometric object of order 2 and it 
is possible to apply the formula which defines the correct variation of energy if we introduce also 
an electromagnetic correction term r (related to the representation chosen) in the definition of the 
variation of the Hamiltonian. This allows to define correctly the energy for the system. Moreover 
we obtain the symplectic structure and the phase space of the theory, with results in full accordance 
with |3|. 

To perform the integration of the variational equation which defines the variation of energy for the 
electromagnetic field we choose two different sets of boundary conditions, corresponding respec- 
tively, to an adiabatic system and to an electrically isolated system [37]. In both these cases we 
obtain a priori a suitable energy contribution to the Einstein-Maxwell Lagragian due to the pure 
electromagnetic field. In the first case this contribution vanishes, this meaning that the contribution 
to the energy comes out only from the pure gravitational part of the Noether charge and it keeps 
track of the electromagnetic field only through the solutions of Einstein-Maxwell field equations. In 
the second case the electromagnetic contribution to energy turns out instead to be proportional to 
the product of the electric charge and the electrostatic potential (integrated over the boundary), a 
result which is physically reliable if compared with Classical Electrodynamics 



Gluing together the results obtained for the gravitational field and for the electromagnetic field, 
we can calculate the energy content of a spatially bounded region of spacetime in the framework of 
Einstein-Maxwell theory and it is then possible to apply the theory developed so far to the case of 
rigidly rotating horizons (see |§, 0). 

The classical definition of black hole thermodynamics, based on the definition of entropy for a 
Killing horizon, deals in fact with quite unphysical models, as it applies only to static or quasi- 
static spacetimes (which means small perturbations from a static situation and thus no radiation is 
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admitted nearby the horizon). Moreover, to define the concept of event horizon for non-stationary 
spacetimes, we need to know the whole history of the spacetime and this is in contrast with the 
concept of physical observer. 

A generalization of these concepts to more physical situations has then been proposed by Ashtekar 
and his coworkers in ||2|, Q (and references therein), where they introduce the notion of isolated 
horizon as a 3-dimensional null hypersurface A embedded into spacetime. Cross sections of isolated 
horizons are, roughly speaking, non-expanding surfaces, isolated from the outside and with a null 
flux of matter and radiation through or outside them. To define the variation of energy and the 
conserved quantities of isolated horizons in the Noether framework we do not need that the spacetime 
has a global Killing vector field but it is sufficient to assume the existence of a local Killing vector field 
for the 3-metric of the horizon; this geometric requirement is fulfilled by isolated horizons. Indeed 
the horizon Killing vector field ensures that isolated horizon are (quasi) locally in equilibrium, but 
they are allowed to admit nearby radiation. 

The Noether formalism developed to define quasilocal conserved quantities for Einstein-Maxwell 
theory naturally applies to isolated horizons and in particular we can define the area Aa, the 
angular momentum Ja and the charge Qa of the horizon through integrals on the cross sections 
At of A and they are conserved on the whole horizon. From the definition of the variation of 
energy we can obtain a first principle for rigidly rotating horizons, which are defined as (weakly) 
isolated horizons with an internal symmetry, generated by a vector field tangent to the cross sections. 
This requirement together with the definition of (weakly) isolated horizons, ensures that when we 
evaluate the variation of energy on a cross section A^ , a first law of black holes thermodynamics is 
defined under the form: 

'^^A = ^'^^A + $(/)5Qa + ^(1)6 Ja (1) 

where A;(;), ^(z), are parameters of the horizon related respectivly to its temperature, its elec- 
trostatic potential and its angular velocity. 

This paper is divided into six Sections. In Section 2 we review the definition of the Hamiltonian 
structure of a field theory and the definition of the variation of energy, via Noether theorem. The 
formalism used is a pure geometric approach, based on the definition of Lagrangian field theories 
on jet-bundles |^6|. In Section 3 we apply the formalism to the case of General Relativity and we 
find an explicit expression for the variation of energy of spatially bounded gravitating systems. We 
introduce here all concepts and notation that are useful for a (3 + 1) approach to field theories. 
In Section 4 we apply our definitions to the electromagnetic field and we analyse the definition 
of variation of energy for Einstein-Maxwell theory. The general theory is applied to calculate the 
energy for two particular sets of boundary conditions, corresponding to different physical situations. 
In Section 5 we finally introduce the boundary and geometric conditions which define an isolated 
and a rigidly rotating horizon. The direct evaluation of the formula defining the variation of energy 
on such surfaces is nothing but the first law of thermodynamics for rigidly rotating horizons. 

2 Geometric framework 

The definition of the Hamiltonian for a field theory, as well as the definition of energy as its on shell 
value, can be based on Noether's theorem. It has been shown in ||2^ that the formula obtained 
via this definition reproduces, in applications, the same results obtained in the ADM Hamiltonian 
formalism ||^, or obtained using a Hamilton- Jacobi analysis of the Trace-K action functional 
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IjToj, [0]. The advantages arising from Noether formalism are related to the fact that the Noether 
(covariantly) conserved quantities are independent on the addition of boundary terms to the La- 
grangian and consequently these latter terms do not influence the definition of energy. In the original 
analysis of Brown and York [^] different boundary terms for the action functional are dictated by 
different boundary conditions and by different choices of the control modes of the boundary data. 
Before entering into the detais of the matter, we shall shortly summarize the geometric framework 
we shall need in the rest of the paper (we address the reader to [22|, |36|, ||3^ for a deeper and more 
rigorous mathematical exposition). 

The configuration bundle for a Lagrangian field theory on a manifold M, with dimension dim(M)=m, 
is a bundle {B,M;7r). A field configuration is a section a of B; i.e. a map a : M ^ B. We can 
choose fibered coordinates {x'^,y^ 



on B and, in these coordinates, the field configuration can be 
locally represented as : ix^,y^ = ct*(x^)). The jet prolongations of the configuration bun- 

dle are denoted by {J^B,M;j^iT). Local fibered coordinates ■■■,ylt^,,^i,) can be chosen on 

j'^B and they denote the spacetime derivatives of fields up to order k. The prolongation of a field 
configuration o" to a jet bundle is denoted by j^a. 

A vertical vector field is a section X of the vertical bundle V{B) = Ker(r7r), namely a vector- 
field which is everywhere tangent to the fibers; locally it can be written in fiber coordinates as 
X = X^-^ = ^y^^ ^''^d it describes the variation 5y* of the dynamical fields. Its prolongation 
to the /c-order jet bundle of B is denoted by j^X and describes the variation of the fields and of 
their derivatives up to order k. If we denote by $t the 1-parameter flow generated by X on i?, its 
prolongation to the jet bundle J^B defines a flow on it, coherently denoted by j^^t- 
The variation of a generic morphism R : J^B — > A"M (where A^M denotes the n-form bundle over 
M for n < m), along the flow of X at any given section a can be defined as: 



t=o dt 



t=o 



(2) 



where dt = $t o o" denotes a 1-parameter family of field configurations on B obtained by dragging 
a along the flow of X. 

A A;-order Lagrangian is a morphism from the A:-order jet bundle J'^B to the bundle A™'(M) of 
volume m- forms on M and locally it can be written as L = C{y''', y\ , ••,?/!,, „ )ds, where £ is scalar 
density called the Lagrangian scalar density and ds = dx Adx A..Adx"^ is the local volume m-form 
on M (we do not assume an explicit dependence of the Lagrangian on spacetime coordinates x^ 
because this dependence is forbidden in natural theories, see 12^ 
The action functional 



(J 



'L 



(3) 



D 



is defined by integrating the pull-back {j^a)*L of the Lagrangian L over a compact region D C M 
with regular boundary dD. According to Hamilton's principle, field equations are obtained by 
imposing the action A_d((t) to be stationary along the flow generated by any compactly supported 
vertical field X € V{B). In this way one obtains Euler-Lagrange field equations (E o j'^^a) = and 
sections a of B which satisfy this variational equation are called critical sections (see [^] , and 

ii)- 



The variation of the Lagrangian along the 1-parameter fiow of the vertical vector field X can be 
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defined as a global bundle morphism 5L : J^B V*{J^B) ® K^{M) by: 

(4) 



< 5L o j V I j^X >=j^{lo o j^a 



t=o 



where V*{-) denotes the dual bundle of the vertical bundle V{-) and < • | • > denotes the natural 
duality between V*{-) and V{-). 

In this paper we consider only natural theories, namely those field theories for which the configu- 
ration bundle B is natural, i.e. each spacetime vector field ^ = £,^d^ lifts naturally (i.e. preserving 
the commutators) to a unique vector field ^ over the configuration bundle as well as it uniquely lifts 
to any prolongation j'^B. This means that for each ^ = ^^d^ G X(M) we have a naturally lifted 
vector field £ X{B) defined by: 

i = e{xnd,+e{x^,y')di (5) 

while j^^ is naturally defined as the prolongation of the vector field ^ to the jet bundle j'^B. 
We stress that the coefficients ^* in (P) depend on and their derivatives up to an algorithmically 
computable finite order r, which depends on the case considered, called the total order of the natural 
bundle. The sum s = k + r of the order of the Lagrangian with the order of the natural bundle is 
called the total order of the theory. 

A Lagrangian theory is natural or, "physically" speaking, covariant if each spacetime vector field ^ 
is an infinitesimal Lagrangian symmetry, which is equivalent to state that the following holds: 

£^L =< 6L I j''£^a > (6) 

where £^a is the Lie derivative of a section a defined by 

£^a = To{i)-ioa = {£i,y')di (7) 

From this formula it is clear that £^a is a vertical vector field over B and it describes the evolution 
of the field configuration a along the fiow of ^. In local coordinates formula (^) is equivalent to 
require that: 

It is well known that General Relativity with the Hilbert Lagrangian is a covariant theory, i.e. it is 
natural; we will show in the sequel that also electromagnetism with the Einstein-Maxwell Lagrangian 
can be treated as a natural theory, provided the configuration bundle and the Lie derivatives of the 
electromagnetic vector potential are defined in an appropriate way |23], [^], |38]. 



The variation bL of the Lagrangian splits as follows by a covariant integration by parts: 

< bL o j^G I j^X >=< E(L) o f^a \X>+d< F(L) o j^^-V | j'^-^X > (9) 



where E and F are called respectively the Euler- Lagrange and the Poincare -Cartan morphism 
We stress that the Poincare -Cartan morphism is not unique for theories of order higher than 2 
(even though in natural theories uniqueness is achieved, for any order k, through the introduction 
of a dynamical spacetime connection; see and [^] for details). When evaluated on a particular 
configuration and on a vertical vector field X, the morphisms E and F are respectively identified 
with a m-form and a (m — l)-form over M. In local coordinates bL can be expressed as: 

bL{j'^a,j''X) = E,{f'^a)X'ds + ¥\f'^-'a,j''-'X)dsx (10) 
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Substituting this expression into the variation dxA£){a) of the action functional (^), it is easy to 
see that Ej represent the Euler-Lagrange equations of motion, while ¥^ is a boundary term which 
vanishes if suitable conditions are imposed on X (usually one states that X is vanishing on the 
boundary dD together with all its derivatives, which is equivalent to state that it is "strongly" 
compactly supported). 

If we deal with natural field theories, Noether's theorem allows us to construct a covariantly con- 
served current and a conserved quantity for each vector field ^ on M. This object can be globally 
constructed and well defined from a geometric point of view ||2^. The Noether current, for a 
solution a is locally defined as a (m — l)-form over M: 

£\L,tcT) = ¥\f^~'a,j^-\£^a)) - eCU'T) (H) 

It is well known that there exists a unique and global decomposition of the Noether current as 
the sum of a (m — l)-form if^ of M (called the reduced current, which vanishes on shell) and the 
divergence of a (m — 2)-form U, named the superpotential of the theory (see p^): 



£{1, a) = £\L, e, a)dsx = £^1, ^, a)dsx + [d^U^f'iL, ^, a)]dsx (12) 

Both the superpotential and the reduced current are linear in the components of the infinitesimal 
simmetry generator ^ and their derivatives up to order (fc + r — 2) and {k + r — 1) respectively, where 
fc + r is the total order of the theory (see our discussion about formula (^)). 

Covariantly conserved quantities are naturally defined as the integral of the current f on a (m — 1)- 
dimensional submanifold T, C M with a compact boundary dT, C S C M: 

Qj:{L,^,a)= [ £\L,^,a)dsx= f £HL,^,a)dsx + l [ U''\L,^,a)ds^x (13) 

When evaluated on shell, since £ is proportional to the equations of motion, the conserved quantities 
are pure boundary terms integrated on 5S: 

Q^{L,i,a) = \f U^''{L,C,a)dsxf, (14) 

which means (see (0)) that the conserved Noether currents are exact on-shell. When evaluated for 
a particular theory and for a fixed ^ these quantities do not reproduce the expected physical values 
because of the "anomalous factor problem" |E^. The solution to this problem relies on the so called 



ADM covariant formalism |£3|, |2^]. According to this, we compute the variation of the conserved 



current (11) along a vertical vector field X as: 

6x£{L,^,a) = 5x¥{L,£^a)-i^{5xL) = 

= 5x¥{L, £^a) - ii:E(L, X) - ii:[dF{L, £^a)] = 

= 6x¥{L, £^a) - £^¥{L, X) - igE(L, X) + d[ii:¥{L, i^^a)] = 

= Lo(L,X,£^a) -ii.¥{L,X) + d{ii:¥{L,X)) (15) 

where 6x is defined by formula (^) and we shortly set ¥{L, £^(7) = ¥^{j'^^~^a, j^~^{£^a))dsx- The 
symplectic current uj (see [jl^) is a (m — l)-form over M 



uj{L, X, £i:a) = 6x¥{L, £^a) - £^¥{L, X) 
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(16) 



which suitably defines a symplectic structure for the field theory. The covariant ADM method 
consists in defining the variation of the conserved quantity, pushing the boundary terms appearing 
in the right hand side of (^) into the definition itself, namely: 

6xQi:{L,tcT) = [ Sx£(.L,^,a)- [ t^¥{L,X)= (17) 

= [ Sx£{L, [ [6xU{L, - ^€F(L, X)] = 

= [ Lo{L,X,{£^a))- [ ii:E{L,X) 
Jt. Jt 

This definition generalises to all natural theories the analysis originally given by Regge and Teit- 
elboim for the ADM Hamiltonian with asymptotically flat solutions and gives the physically 
expected results for all the conserved quantities of the theory. 



The Hamiltonian structure of the theory naturally arises from this definition P], |23], [28|, |37]. 



We define the variation of the Hamiltonian as the variation of the conserved quantity, relative 
to a Cauchy surface S. The fundamental requirement is that the infinitesimal generator ^ of the 
symmetry, which defines 6H, is transverse to the surface S. From equation (^^ we have that: 

dxH{L,^,a) = [sx£{L,ta)-[ i^¥{L,X) = 
Jt Jot 

= [ uj{L,X,£^a)- [ i^'E{L,X) (18) 
Jt Jt 

The variation of the energy is defined as the on-shell variation of the Hamiltonian Q, [37|; in this 
case E(L, X) vanishes and the variation of energy is defined as the integral of the symplectic current 
uj on S. We remark that if both X and £^a are solutions of the linearized field equations, from 
([l6|) we have that u! is closed on-shell ]l3|, [|l9|] . 

The Poincare -Cartan form is linear in X and its derivatives up to order {k — 1) and is linear in 
and its derivatives up to order (k + r — 1), where r is the order of the natural bundle. When we 
integrate (^) on S it is then possible to split the symplectic current into two terms (through an 
integration by parts): 

/ uj{L,X,£^a)= [ cD(L,X,i:^a)+ / d[T{L, X, £^a)] (19) 
Jt Jt Jt 

where a) is a (m — l)-form, while r is a (m — 2)-form which can be integrated on the boundary dTi 
using Stokes' theorem^. Our aim is thence to suitably correct the definition of the Hamiltonian and 
consequently of energy in a way that the symplectic structure of the theory is completely defined 
by Jj^oj, i.e. by an integral over the (m — l)-surface S, while energy is a pure boundary term, i.e. 
an integral over dT,. 

^ It is the analogy with Classical Mechanics which suggests how to perform the splitting: roughly speak- 
ing, our purpose is to generalize to natural field theories the well-known formula SH — qSp — p5q for Hamil- 
ton equations in Classical Mechanics. Hence, in the integral f^u>{L,X, £^a) all the terms which are not of the 
form Jj, [{£^(j) Sxp — {£^p) Sxo'] d^x have to be pushed, through integration by parts, into the boundary term 

jg^T{L,X, £50-). 
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According with the prescription of the covariant ADM method we can redefine a corrected variation 
of the Hamiltonian 5xH by pushing the new boundary term r into the definition of 6xH itself, i.e.: 



6xH{L,C,^) = 6xH{L,^,^)- [ T{L,X,£^a)= (20) 

JdT, 

= [ 6x£{L,^,a)+ [ [5xU{L,0-^¥{L,X)-T{L,X,£^a)] (21) 

The new variation of the Hamiltonian does no longer contain any boundary term. In fact boundary 
terms arising from the variation of the reduced current 6xS are completely cancelled by the variation 
of the superpotential together with the correction terms and r, so that, from equations ([l8|) 
and (p^), 5xH finally results to be: 

SxH{L,^,^)= [ Cj{L,X,£^a)- [ ^E{L,X) (22) 

We can now define the variation of energy as the on-shell value of the variation of the Hamiltonian 



. It turns out to be: 



6xE{L,i)= [ [6xU{L,0-HnL,X)-T{L,X,£^a)] (23) 

because 6x£ in ([21] ) is identically zero since we have assumed that X be a solution of the linearized 
equations of motion. This master formula ( p3D gives us a recipe to define the energy E(L, S,) once 
the variational equation dxE is solved. We stress that the definition (|2l] ) of the variation of the 
Hamiltonian and consequently the definition (^3|) of energy variation 5E do not depend on any 
divergence term possibly added to the Lagrangian. 

We indeed remark that the final formula of 6H is invariant under a change of the Lagrangian by 
pure divergence terms. This means that the variation of the Hamiltonian and the variation of energy 
are defined for an equivalence class [L] of Lagrangians, rather than for a single one, defined by the 
equivalence relation: 

L' ^3T £ A™-^(M) =^ L = L' + dT 

which implies that L and L' generate the same field equations. Alternative definitions of energy 
that can be found in literature define 5H by adding some boundary terms to the Lagrangian or to 
the canonical Hamiltonian ||^], |31]. These boundary terms which do not affect the equations 
of motion modify the definition of energy and consequently the control mode of the fields on the 



boundary [37|. Suitable terms have then to be added a posteriori to the action functional to obtain 
the predefined boundary control mode. 

The symplectic structure of the theory and consequently its phase space are instead defined in 
literature starting from the symplectic current ui defined on the surface S and by a symplectic 
structure on the boundary 3S, which is dictated by the choice of the boundary term we add to the 
Lagrangian; see [^. 

The definition of the variation of the Hamiltonian we propose gets rid of all these problems, since the 
symplectic structure is defined only by the surface part Jj, uj of the symplectic current. In analogy 
with Classical Mechanics, as shown in [^, |28|, [37|, the variation of the Hamiltonian contains a 
term which determines the symplectic structure of the theory as well as a term which contains the 



Lagrangian equations of motion. This is just the case of (p^), where a) is related to the symplectic 
structure while the second term is linear in the Lagrangian equations of motion so that, eventually, 
formula (21) can be considered as a well-defined expression for the variation of the Hamiltonian. 
Moreover, boundary conditions in this framework do not play any role either in the definition of 
the variation of energy nor in the definition of the symplectic structure of the theory. They will 
instead play a key role in the integration of the variational equation 5xE where imposing boundary 
conditions means to fix the values of the vertical vector field X (i.e. of the variations of fields) at 
the boundary dT,. The choice of boundary conditions, namely the choice of the vertical vector field 
X, together with a reference "zero level" for energy, provides us with the definition of energy |1C], 

m- 

We will see that this formalism applies to Einstein gravitational field and also to the electromagnetic 
field in the framework of a geometric natural formulation of Einstein-Maxwell theory. This is possible 
as the total order of the electromagnetic theory results to be s = 3. 

The theory has been formulated to calculate energy for a bounded, compact region of spacetime. 
However it is possible to consider the case of a system of infinite spatial extent, i.e. dT, in ( p3| ) 
is now the spatial infinity. Results obtained using the quasilocal definition of conserved quantities 
reproduce the results obtained through the standard ADM definitions. Moreover the definition of 
the variation of energy (^) is homologically invariant, if ^ is a global Killing vector field |25]. In the 
case that spatial infinity is homologic to a finite bounded surface B, conserved quantities evaluated 
on B correspond to the total conserved quantities at spatial infinity. 



3 General Relativity 

The main application of the formalism exposed in the previous chapter is General Relativity in 
vacuum with dimM=4. The Lagrangian of the theory is the Hilbert Lagrangian: 

Lh = ^Vag'^'Rf^uds (24) 

where fc = Svr in geometric units, with G = c = 1. This argument has been analysed in a previous 
paper ||2^ (for details on calculations and notations see also [|lO|); we give here just a brief overview 
to recall the main results. We consider a four dimensional manifold D C M of spacetime and we 
assume D to be diffeomorphic to the product S x J?, here S is a 3-dimensional closed manifold with 
boundary 5S = B. For any t & M there exists an immersion: 

ft-^^^t (25) 

and we require T,t to be a portion of a spacelike Cauchy hypersurface. The set of all leaves defines 
a foliation of D in spacelike hypersurfaces. Each St intersects dD in a two dimensional surface Bt 
which is diffeomorphic to B. We denote by B the timelike hypersurface B = dD = IJt -^i- 
The time evolution is generated by a vector field ^ transverse to S^; we can require that CV^t = 1 
and ^ is tangent to the boundary B. We denote by the timelike future directed normal to St and 
by the outward pointing spacelike normal to Bt in S^. The vector field ^ can be decomposed as 

= Nuf" + Nf" (26) 

where the shift lies in TS^, i.e. N^u^ = (this matter will be discussed later in detail). 
We define n'^ the outward pointing spacelike unit normal to B while by we define the timelike 
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future directed normal to Bf in B. Barred and unbarred vectors are related by boost relations (see 



my- 



where we have set: 



jv = Ufj,n^ = —n^u^ = sinh{9) 

V = = 

N 



(27) 



The parameter 9 is usually called the velocity parameter, while v is the boost velocity and 7 = 
(1-t;2)-V2. In fact, if we consider a point on Bt at a time t we can observe its evolution according 
to the unboosted (unbarred) observers, which means observers at rest on Sf and evolving with 
the normal vector u^. Otherwise we can describe its evolution with respect to a boosted (barred) 
observer evolving with the four velocity and seeing the vectors u'^,n^ as normals to Bt. The 
scalar v is in relation to the boost radial velocity between the two classes of observers, while 9 



describes the non-orthogonality of the foliation; see |1C] 



The metrics induced on the surface and B by the projector operators are defined respectively as: 

while the metric on Bt can be defined with respect to boosted or unboosted observers as: 

(^iiu = 9iiv - riiinu + Ufj,Uu = gfiu - n^n^ + u^u^ (28) 

In the sequel we denote by g,h,^,a the absolute values of the corresponding metric determinants. 
The extrinsic curvatures K^i, of St in M, Q^^, of i3 in M and /C^jy of Bt in Sj are defined respectivly 
as: 

K — — ft"V V 
KLfip — c ^Df^rii/ 

where D is the metric covariant derivative with respect to the metric /i of Sf . We denote by K, 
0, /C the above quantities contracted with the corresponding metrics. We shall instead denote by 
ICf^i, the extrinsic curvature of Bt with respect to the normal n'^, i.e. }C^,y = ^K.^^ — Iva'^OuV aUp] 
see [p^] . 

The three momenta P^'^ of the surfaces T,t and H'^'^ of B are defined in terms of the extrinsic 
curvatures of the 3-hypersurfaces in M: 

The time evolution vector field can be decomposed, into a normal part to called lapse and a 
tangent part N^^ called shift, by: 
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Equivalently it can be decomposed with respect to the boosted observers: 

where N = ^ and N^' is the projection of N^' on Bt, i.e. = a^N°'. 
3.1 Variation of the Hamiltonian in General Relativity 

It is now possible to calculate explicitly the symplectic current and the variation of the Hamilto- 
nian for General Relativity in vacuum, which is governed by the Hilbert Lagrangian (details and 
calculations can be found in [p^]). 

On a leaf Sj of the foliation, the reduced symplectic current uj and the boundary part r in (|l9|), 
turns out to be: 

Co{LH,X,£^g) = [{£^h^,)5xP^'' -{£iPn^xKu]d'x (29) 
T{LH,X,£^g) = —[£^{^/an^5xu^) - 5x{\/^n^£^u^)](fx (30) 

From these expressions for the reduced symplectic current Co and r it is clear that the symplectic 
structure and the phase space of the theory are both completely determined on the Cauchy surface 
Tit in terms of the 3- metric h^^ and its conjugated momentum P^'^. 
The variation of the Hamiltonian ([2l| ) has been calculated in [^] and turns out to be: 



6xH{Lh, e, ^t) = [ {HbxN + ?^„<5xiV" + [h^,\6xP'"' - [P^''\6xKu}(f'x (31) 
where [/i^j;/]^ and [P^'^^ follow directly from the 3-dimensional Einstein equations. Notice that no 



boundary term appears in (|3l|). Comparing (l29| ) with (31) the Hamiltonian equations (|2^ ) read as 
follows: 

-^G^'^u^u, = n = Q 

^ ^G"^/i^/i^ = £^P^"' - [P^"']^ = 
For our later purposes it is not necessary to formally integrate (|3l|) and to give the explicit expression 



of H (see |10] and |28| for details). Using the formalism developed, in fact, one can directly turn 
to calculate the time rate of change of the Hamiltonian. It is easy to show that the Hamiltonian is 
conserved along the flow of a vector field iff ^ is a Killing vector for the boundary metric on B, 
i.e. £^j^u = 0. 

The variation of the energy turns out to be equal to (see again |2^] ): 

6xE{Lh,0 = [ {Ndx{V^e)-N''6x{V^fa) + ^^s''''Sxa^,A(fx+ (32) 

J Bt 

+ y[ [£^{^)5x{9)-6x{V^)£de)]cfx 
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where we have set: 



-e=m 



To integrate the expression ( p^ ) we need to impose suitable boundary conditions. C.-M. Chen 
and J. M. Nester have shown that there exist only two different Hamiltonian boundary terms 
which correspond to covariant boundary conditions and they are respectively the Dirichlet and the 



Neumann control, which define two different energies |14|. The energy corresponding to Dirichlet 
boundary conditions is usually accepted as the internal energy; Dirichlet boundary conditions fix 
the 3-metric 7 on the boundary B: 

5xN \b= SxNf" \b= Sxct^^u \b= (33) 

which can be considered as a restriction on the vector field X and its flow. This last conditions 
allow to integrate (|3^) and we obtain: 



E{LH,g,Bt)-Eo{LH,go,Bt) = 1^ [V^{Ne - N^j^) + ^£^{V^){e)}(fx (34) 

We stress that we are integrating along a 1-parameter family of solutions all satisfying the same 
fixed boundary conditions (|33[). The term Eo{gQ,Bt) corresponds to the quasilocal energy of a 
background solution go inside the 1-parameter family and it becomes the "zero level" for energy. 
Owing to the properties N = Nq |g, iV" = Nq |g and a^i, = (jQ^y |g, we obtain the equivalent 
expression: 

E{LH,g,Bt) = {V^[iV(e-- eo) - iV'^(ja - /oa)] + \£d^)iG " ^o)]d\ (35) 

If we assume the vector field ^ to be a Killing vector on the boundary, this guarantees the conserva- 
tion of energy in time and the term £^{y/a) identically vanishes, so that formula ( |35|) agrees with 
the definition of energy given in [10|. 



4 Electromagnetic field 

In this Section we shall describe the Einstein-Maxwell theory and we shall calculate its energy in 
relation with different boundary conditions. The energy of the system is related to the control-mode 
we choose. It is possible to imagine a "gedanken experiment" which reproduces in a laboratory the 
same conditions described; controlled quantities and energy have in this case an interpretation as 
physical osservables , [^] . 

The approach we use here is different from the the ones usually found in literature. We treat 



electromagnetism, described by the Maxwell Lagrangian, as a natural theory (|23] and [^), which 
implies that we consider it as a [/(l)-gauge theory based on a natural principal bundle. This choice 
implies that is possible to construct a fiber bundle of geometric objects on spacetime (describing 
the configurations of the electromagnetic field) which is a principal bundle over M. An immediate 
consequence of this fact is that the magnetic charge of the configurations considered is identically 
zero. This is the same case usually treated in literature to describe the Hamiltonian formalism 
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for Einstein-Maxwell theory using a trace-K action but in this latter case the vanishing of 

magnetic charge is assumed (somewhat equivalently) as a restriction on the class of solutions. 



4.1 Einstein-Maxwell theory 

We consider a Maxwell theory described as a gauge theory with gauge group U{1), based on a 



fibred principal natural bundle (P, M, U{1);tt), whose existence has been proven in p4|| and |25]. A 



configuration is a section of the connection bundle ^/^jy ) which we denote hy A : M ^ 'ulTJ' ^^^^^ 
A = A^{x)dx^; it is called the quadripotential of the theory. The coefficients A^ are M(l)-valued, 
where m(1) = iM is the Lie algebra of the Lie group. This means that can be interpreted 
as a ?7(l)-connection and its curvature F^i, is the field strength. We are particulary interested 
in the interaction between electromagnetic fields and gravity, with configuration bundle ( '^u{i) ^ 
Lor(M),M;7r) and described by the Einstein-Maxwell Lagrangian density: 

J^E&dfgJ^A) = Cnifg) + ^MigJ^A) (36) 

with: 

CMigj'A) = -l^^g-^^gP-F^pF^, (37) 
We have defined -F^^ as the curvature of the quadripotential 

F = dA = A dx^ = ]^{dcAp - dpAc,)dx'^ A dx^ (38) 

while the naive momentum with respect to the curvature F is denoted by: 

= -k^ = Vgg'^^g^'F.u (39) 

The components of the curvature F can be decomposed on the Cauchy surface T,t, hence defining 
the electric fiux density and the magnetic induction of the electromagnetic field: they represent the 
physically observable quantities. The Lie derivative of the electromagnetic potential A^, according 



to ||2^ and |25|, is defined as: 



£^A^ = C{d,A^) + A,{d^e) + qd^i^n = A^ + qd^{d,n (40) 

where A^ is a shortcut for A^ = ^^[d^A^) + A^^d^S,^) and g is a dimensional constant which 
represents the geometric charge of the field carried by the natural representation 

Rq-.Xe GL{m,R) exp(^z^ | det{X) \ ) e U{1) (41) 

(here e < is the unit charge of the electron). This representation (^ij) defines the principal bundle 
(P, M, ?7(1); tt) as a bundle associated to the frame bundle (or in other words it is the representation 
which defines (P,M, [/(!)) as a natural bundle; see |2^, [^]). Notice from (^) that the bundle 
of connections TT^Ty turns out to be a natural bundle of order r = 2. Indeed, the Lie derivatives 
£^A of its sections depend on the coefficients .^^ and their derivatives up to order two. Being the 
Maxwell Lagrangian (^^ of order k = 1 the total order of the electromagnetic theory turns out to 
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be s = 3. 

The Maxwell Lagrangian will be considered separately for the sake of convenience; the generalization 
to the Einstein-Maxwell theory can be easily obtained recalling the results obtained in Section 3. 
Euler-Lagrange equations of motion follow for Cm immediately from the variational principle ( |lO|) : 
we obtain that the electromagnetic field must satisfy the equation J'^ = 0, where 

J'^ = -^(9„r°) = -^(V«r°) (42) 

is called the current density. We remark that the Poincare -Cartan form is obtained by the pre- 
scription (P) as: 

¥\LM,X) = jf>'HxA^ (43) 
k 

Using the definition of the Noether current, from formula ([TT|) it follows that: 

£\LM,0 = \f^^{£^A^)-eC (44) 

and from ( |40| ) we obtain the natural splitting of the current into a reduced current term and the 
divergence of the superpotential 

£^{LM,i) = £\Lm, + d^U^^'{LM,0 (45) 
where in local coordinates we have 

having defined the stress-energy tensor of the electromagnetic field by setting: 

= -Ig'^iF^F^, - \g^.F^PF^p) (46) 

The Lagrangian of the Einstein-Maxwell theory is the sum of Lh and Lm = Cuds. Accordingly 
Noether currents are simply obtained for this theory as the sum of ( ^5| ) with the analogous expression 
for General Relativity (see 



4.2 Energy for the Electromagnetic Field in a (3 + 1) perspective 

We are now going to apply the definition ( |23| ) for the variation of energy to obtain an expression 
for the Einstein-Maxwell Lagrangian, writing for convenience only the terms regarding the electro- 
magnetic field, since the ones for the gravitational field are well known from the previous Section. 
To obtain the final result we need to integrate the superpotential and the correction terms r on 
the boundary Bt of the spacetime slice S^, which in turn is assumed to be a Cauchy surface. We 
shall consider a (3 + 1) splitting of the dynamical fields evolving with respect to boosted observers, 
which means that we put ourself in a reference frame which corresponds to observers evolving with 
respect to ^. 

We define the electric vector field E and magnetic vector field B in the usual way: 
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where ei^^^^ is the totally skew-symmetric Levi-Civita tensor on M. Conversely the curvature can 
be defined in terms of the electric and the magnetic vector fields as: 

Fa(3 = {Ec^up - Epua) + u^e^afif.B'' (47) 

We define the electrostatic potential and the normal component of the electric field as: 

$ = -Aau" 

Analogously, by {A^,E^,B^) we denote the quantities projected on the 2-dimensional surface Bt, 
so that, e.g.: 



(48) 



(and analogous expressions for the other quantities). The symplectic current defined in (|l9|) owing 
to (^) turns out to be: 



(49) 



This expression justifies the definition of f^^ as the conjugate momentum with respect to the 
electromagnetic field. The pullback of the symplectic current (^) on the surface can be written 
as: 



(50) 



u:{LM,X,£^a) = [5xS^{£^,A^) - £^S^{5xA^)]d?x = 
= [5xS^[£t,A^)-t^'{5xA^)]d?x 

where we have defined the 'i- dimensional electric density £^ and its time rate of change by: 

= IVhEf" 

Using the formula (40) and recalling that the total order of Maxwell-Einstein theory is s = 3, we 
can split /j^w, in accordance with (p^), into two bulk terms and a boundary term: 



St 



[5x£nA^)-£^{5xA^)y^x + 



(51) 



q\ / 5x{D^£^')d,ed''x + 6x 



Finally it is possible to identify the reduced symplectic current and the term r (see (|T9D): 

Cj{LM,X,£^a) = [6x£''{A^)-£^{dxA^)]d^x 



t{Lm,X, £^(t) 



2q 



6x{E^n^V^}{dpe)d'x 



(52) 
(53) 



while the remaining term in (51) is identically vanishing on shell. This last expression encodes the 
symplectic structure of the theory and defines both the phase space and the conjugate momenta 
[13|. We remark that the boundary part t{Lm,X, £^a) of the symplectic current is non vanishing 
since this holds whenever the total order of the theory is higher than 2: for electromagnetic theory, 
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treated as a natural theory, the total order is in fact s = 3. This term tm will influence the definition 

of the variation of energy in accordance with (23). 

The reduced current ( [45| ) for this theory can be rewritten as: 

£{Lm,0= - {iNTJ'uau'') + {T^uxN'')}Vhcfx + [gApCP]Vhd^x+ (54) 
+ g[n,>5^)a,r-a,(j["<5^1)r^^a]\/^rf^a; (55) 

where we have set: 

g = -DpS^ (56) 

which is the spatial part of Gauss' constraint and vanishes on-shell. If we define: 

(n^ = -Vh{T,^uxun = ^{E^Eo^ + B^B^) 
\ n^o. = Vhh^aiT^^ux) = ^{e^xpE^B") 

we finally have that: 

/ iiLM,0= f [Nin^ -^g) + N''in^^ + A^g)]d^x+ (57) 

+q [ Kj("5^)a,r-9.(j[°5^i)r«jv^d'x 

To calculate the variation of the energy by means of (]23| ) we calculate now the expression for the 
superpotential, evaluated on Bf: 

[ U{LM,Od^x = r f E^[-N<^> + ApNP + vN{A^n'')]^/^(fx+ (58) 
JBt k Jb, 

- q6x [ J^''5SunaC^d^x + ^6x [ {Ef^n^dp^^jd^x 
together with the correction term, which turns out to be: 

/ ^¥{Lm,X) = j[ NiE^u^" + vE^" + usn^e^''''^Bp)5xAp^d^x (59) 

JBt JBt 



Neglecting terms which vanish on shell, from (|5^), ( pq ) and (pOD we obtain that the variation of 
energy for the electromagnetic part of the Einstein-Maxwell theory is: 

5xEiLM,C) = [ [5xU-^¥-T]d^x= (60) 

JBt 

= -\ [ N[^-v{n"A^)]5x{E^y/^)d^x+'^ [ N"5x{E^ A^./^)d^x + 

k JBt ^ JBt 



Bt 



Bt JBt 

iV[u^" + Uynpe''^''^B^]5xAo,y/ad^x 



We will see that the first term in the last expression is directly related with the charge of the field, 
the second is related with the rotational degrees of freedom and the third one is determined by the 
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values of the electric and the magnatic field projected onto the surface B^. 

We remark that this expression of energy exactly cancels the boundary terms in the variation of the 
reduced current ( ]57|) so that, in accordance with the general theory, we obtain that the variation of 
the Hamiltonian (pi]) 

6xH{Lm,^, S) = / 6x£{LM,t ^)+ [ [6xK{Lm,0 " ^cF(Lm, X) - t{Lm,X, £^a)] (61) 

Jt. JdT. 

is a pure bulk term, which, according to ( p2[ ) and ([5^), generates the Hamiltonian equations of 
motion. 



4.3 Energy and boundary conditions 

In the Einstein-Maxwell theory we have that Hilbert Lagrangian is minimally coupled with the 
Maxwell Lagrangian, so that the conserved quantities are additive: 

5E{Lem,0 = ^E{Lh,0 + SE{Lm,0 (62) 
We have previously defined in (35) the quasilocal internal energy for the gravitational field contained 



in a region of spacetime bounded by a surface B, imposing the Dirichlet metric boundary conditions 
(p3|). In analogy, we want to derive now the energy contribution due to the pure electromagnetic 
part of the Einstein-Maxwell theory. To perform calculations we rewrite equation ( |60|) in terms of 
boosted observers: 

5xEiLM,0 = -T [ N^6x{E^V^)d^x+l [ N''6xiE^AaV^)d'^x (63) 
- J I N[u^npe^P^^B^5xA^\^<f 



where we recall from from the definitions given above that E-^ = E-^, $ = — A^u", B = vyeu^apF'^^ 
and Uynpe''l^°'^' = e^^^^n^n^j. 

To integrate this variational equation and consequently to obtain the energy contribution in the 
region bounded by Bt we have to consider a 1-parameter family of solutions of field equations which 
admit the same physical data fixed on the boundary. We can choose different control modes which 
correspond to two different physical situations, namely a system which is electrically isolated from 
outside or an adiabatic system (see p7|]). 



• We can choose a control mode for the boundary components of the electromagnetic potential 
Aa (i.e. we control the magnetic flux through Bt) and for the electrostatic potential which 
is equivalent to state that the laboratory is electrically isolated from outside. In other words 
we assume 6xAa \Bt= ^ \Bt= ^ (|63|). Recalling that Dirichlet conditions ( |33[ ) are 

chosen for the gravitational fleld 6x"f \b= 0, we obtain from (|6^) that: 

6xEiLM,0 = ^^{-1 J N^E^^dPx + {N"E^A^^)cfx} (64) 
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This expression is integrable so that: 



E{Lm,0-Eo{Lm,0 = -t [ {N^-N''A^)E^^d^x= (65) 

k JBt 

JBt 

We recall that the variation is performed along a one parameter family of solutions, all ad- 
mitting the same boundary conditions. Hence Eq is the energy corresponding to a reference 
solution inside this family. 



Expression (65) for the energy of the electromagnetic field, owing the boundary conditions for 
the gravitational and the electromagnetic field (l> |bj= Ib* j ^Oo \Bt = Ist, I/mu = iQ^u)-: 
can be rewritten in an analogous way as: 

^(Lm,0 = -T I (^^ - iV%)(E^ - E^)^dFx (66) 

JBt 

where the subscript clearly refers to the reference solution. This formula reproduces the 
energy content of the electromagnetic field and it stresses that the response variable, with this 
boundary condition, is the normal component of the electric field. The one parameter curve 
of solutions is parametrized by E-^, i.e. by its electric charge. 

• If we choose instead a control mode for Aa and (i.e. 6xAa |_Bt= and 6xE-^ \Bt = 0) we 
state that there is no flux of magnetic and electric field through Bt- In analogy with classical 
thermodynamics this system can be defined to be adiabatic. In this case calculations are even 
more trivial. Prom equation ( |63|) and always assuming Dirichlet boundary conditions (|3^ we 
find: 

6E{Lm,0 = ^ E{Lm,0 = const (67) 

This result, which appears to be quite shocking, can be interpreted by saying that the con- 
tribution to energy coming from the pure electromagnetic Lagrangian is neglectable. This is 
in accordance with the fact that the system is adiabatic. Nevertheless the pure gravitational 
contribution (|35| ) to the total energy keeps track of the electromagnetic field through the 
metric solution g of Einstein equations G^u = fc7^,y, where 7^,^ is the electromagnetic stress 
tensor defined in (^6|). In other words, the quasilocal energy (^) of the gravitational field 
takes account of the presence of the electromagnetic field. 

These results are in accordance with the analysis on control modes proposed by Kijowski in |37], 
using a symplectic type analysis in the framework of Legendre transformation and Hamiltonian 
formalism. We remark that formula (|6^) generalizes this analysis to the case of Einstein-Maxwell 
theory in a spacetime region with non-orthogonal boundaries. 



5 Rigidly rotating horizon 

In this Section we are going to recall the definitions of isolated horizon, weakly isolated horizon and 
rigidly rotating horizon. The concept of isolated horizon was first introduced by A. Ashtekar and 
coworkers in Q and later developed in ^ to obtain a viable and more general definition. This 
new concept allows to generalize the laws of black holes thermodynamics to the case of black holes 
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which do not admit an event horizon and a global Killing vector field. This situation is much more 
realistic with respect to the quasi-static laws of black holes mechanics based on the existence of 
an event horizon and referred to equilibrium situations and small perturbations out of them. The 
definition of isolated horizons is given intrinsecally and this implies that we do not need to know 
the whole hystory of spacetime to define this geometric surface (as we usually need for the event 
horizon in non stationary spacetimes). We can admit isolated horizons with radiation infinitesimally 
near the horizon surface, meaning non-stationary spacetimes. Conditions are imposed only on the 
internal geometry of the horizon, not on the whole spacetime, to ensure that the black hole itself is 
"isolated". This is reasonable also in analogy with the classical approach to thermodynamics ||l5[| 
where equilibrium is just defined in function of the internal degrees of freedom of the system. 
Physical examples of isolated horizons can be found in collapsing stars or in cosmological horizons 
in de-Sitter spacetime (see |^ for a review on the matter); in the latter case no singularities are 
present in spacetime, but it is still possible to define thermodynamics for such surfaces. 
The matter covered in this subsequent Section is (apart from minor changes needed to adapt to our 
notations) a brief account of Ashtekar's definitions of |^], which we consider worth of being recalled 
for the sake of completeness and clarity. 

5.1 Geometry of horizons 

The fundamental requirement is to ask the isolated horizon A (a full definition will be given later) 
to be a null non expanding 3-surface, which means that the 3-metric on it admits a Killing vector; 
this encloses the fundamental concept that A is isolated in a suitable physical sense. If any matter 
or electromagnetic field is in interaction with the gravitational field we have to state that the fiux of 
radiation and matter through the horizon is zero, which for the electromagnetic field is equivalent 
to say that the relevant component of the Poynting vector is zero on A. With more precision we 
can start with the following definition: 

Definition 1 A null hypersurface A in a 4- dimensional spacetime M is a 2>- dimensional subman- 
ifold of the spacetime admitting a null normal vector I = l^d^ (which for definition is also tangent 
to the same surface A). 

In the sequel we will consider only null surfaces which are diffeomorphic to the product A ^ Mx S'^ 
such that for each t £ R there exists a diffeomorphism ipt : S"^ ^ At, where A^ is a two dimensional 
spacelike surface. 

The metric defined on A is clearly degenerate for definition (see p6|). We now consider a region 
D of spacetime such that A C D; like in the previous Section, we assume D to be foliated by 
3-dimensional spacelike hypersurfaces T,t, which induce a preferred foliation ip of A and we denote 
again by the future directed causal normal to Sf and by n'* the outward pointing normal to A^ 
in St. 

To define the metric on A we have to introduce a null vector field / transverse to A. We ask this 
vector field to be normalized relatively to / so that I"' la = —1. Notice that I" £ TA, 1^ G T*A, while 
transversality requires 1^ ^ T*A and P ^ TA. Due to this normalization condition it is possible to 
express / and / in terms of and as: 

la = fiUa + ria) 
la — 2/^^*^ ^a) 
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where / is a non- vanishing function on M. Now the "metric" and the projector operator 
onto A can be defined as: 

QaP = gal3 + 2Z(q,/;3) 



1/3 — 9 (lap 

Since q^p is degenerate the " inverse" g"^ of this " metric" is not unique but any reasonable definition 
of it remains unchanged under the addition of terms of the form l^°'V^\ with V tangent to A. 
The non-degenerate metric a over Af can be obtained by considering either A^ as a submanifold of 
A or a submanifold of Ej. We have: 

I = 5^ + + hp 

where the analytical expression of cJq,^ is the same of qa/3, but has not to be confused with q^/s which 
is defined in a different space. The time evolution, defined by a vector field transverse to S^, see 
(|26|), can be expressed as an element of TA as: 

= ivr + N" (68) 

where iV = ^ and iV" = a^N^^. 

For each null hypersurface A it is possible to select (non uniquely) a vector (which is for definition 
a null normal vector) and define the tensor expansion of A as: 

0(i)^f, = qlq^^yj^. (69) 

The expansion plays for null surfaces the same role that extrinsic curvature plays for spacelike or 
timelike hypersurfaces. As /° is normal to A we have immediately, from Frobenius' theorem, that 
the twist of the expansion is zero = 0, so that O^p can be decomposed as: 

= 2^{l)'^af3 + Ca/3 (70) 

where = C{ap) is called the shear of the expansion tensor. The scalar: 

1 



^w = %).,^"" = 7^^K^/^) (71) 

is properly called expansion of the null surface (see ) . We are now ready to give the definition 
§: 

Definition 2 A non-expanding horizon A C M is a 2>- dimensional submanifold of spacetime {M,g) 
with the following properties: 

• A is a null hypersurface and it is topologically the product of a two-dimensional sphere and 
the real line: A = S'^ x M 



The expansion 6(^i-j of any null normal to the surface vanishes on A 



• Equations of motion hold at A and the stress energy tensor for matter interacting with gravity 
is such that (— Tq^Z") | a is future directed and causal for each future directed null normal I to 
A 
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We remark that the fundamental property to identify an isolated horizon is to impose = 0. The 
topological condition is equivalent to the requirement that there exists a foliation of the hypersurface 
A in "time" constant surfaces. The matter condition is very weak and it is implied for example 
from the dominant energy condition; it is equivalent to state that there is no flux of matter through 
the horizon; see Q for a detailed presentation. 

Since the twist of vanishes, from the Raychaudhuri equation [^2|, we obtain that every null 
normal I is free of expansion, twist and shear on the non-expanding horizon, i.e. |a= 0, 

6(1) |a= and |a= 0. 

The same Raychaudhuri equation ensures that there are some restrictions on the Ricci tensor and 
in particular that the following should hold 

Raf^n^ \a= 

The natural connection 1-form lo is defined on A (since / is expansion and twist free) as: 

qi^V^f = (loJ^) (72) 
which implies the property for / to be a Killing vector of the 3-metric q: 

£lQa(3 = 

The acceleration k(i-^ of I, defined as: 

l^^V^f = k(i/ (73) 

is expressed in terms of the natural connection as k{i^ = Q. Notice that k{i^ is the counterpart 
for isolated horizons of the surface gravity for Killing horizons |46]. 

Finally we stress that from the third condition in Definition 2 it follows that the stress energy 
tensor contracted twice with is identically zero, i.e. T^^/"/^ = 0. This implies that in the case 
of Einstein-Maxwell theory we are dealing with, = F^pl'^cr^ = and = F^pl'^cr^ = 0. This is 
equivalent to state that the electric and the magnetic components tangential to At are identically 
zero. 

Definition 3 A non-expanding horizon becomes a weakly isolated horizon (A, [/]) if it is equipped 
with an equivalence class [I] of null normals satisfying the equation: 

£iLu = Oyie[l] (74) 

where the equivalence relation is defined as: I ^ I' <;=^ / = cl' ; c £ M is a non zero constant. 

We remark that a Killing horizon is automatically a weakly isolated horizon. As it happens for 
Killing horizons it is possible to define (apart from a constant factor, which depends on the rep- 
resentative inside [I]) the surface gravity for weakly isolated horizons if we equip them (as it has 
been done in the definition) with a preferred family of null normals. The value of k^i-^ depends on 
the element of the class we choose. One can otherwise uniquely define /cj-;) by demanding it to be a 
suitably defined function on the horizon; see Q. 

The properties of weakly isolated horizons ensure that the zero law holds true. This is equivalent 
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to say that the surface gravity is constant on a weakly isolated horizon. In fact from the property 
£lUJ = 0, we obtain that: 

d{k(i))W=Q (75) 

We will consider for applications only the case of rigidly rotating horizons interacting with electro- 
magnetic fields, i.e. the Einstein-Maxwell described above. The theory developed hereafter is also 
true in the case of isolated horizons, which can be considered as particular cases of rigidly rotating 
horizons (roughly speaking with vanishing angular momentum). |^ 

Definition 5 A weakly isolated horizon (A, [/]) is said to he a rigidly rotating horizon (A, [/],(/?) if 
it admits a rotational simmetry ip ( with (p tangent to the surfaces A ) with closed and circular orbits, 
such that £,pl°^ = 0, £ipiOa = and £,pqa/3 = 0. 

This last definition is a restriction on the class of weakly isolated horizons we are considering and 
it is analogous to state that the 3-dimensional geometry of the hypersurface is axisymmetric, with 
if as its infinitesimal rotational symmetry generator ||5|. 

In this case, it is possible to choose a foliation of A in a preferred way so that the surfaces 
intersect A at the preferred 2-surfaces Aj topologically diffeomorphic to 5^ and we state that there 
exists a tangent vector ip" £ TAt, which generates the desired symmetry. Correspondingly, the 
time flow transverse to the surface is adapted to A if = T" = P — Q(^i-j(p" on A. 
In the case of interaction with matter, and particulary in the case of Einstein-Maxwell theory we 
are going to analyse, it is natural to require 99 to be a symmetry also for the electromagnetic field 
projected on A, so that £ .^liq^qu Fafj) = 0. 

In stationary and asymptotically flat spacetimes one usually sets the electrostatic potential to be 
^ = —^f^Af^ and requires a gauge flxing on the electromagnetic potential such that tends to 
zero at inflnity and £^Af^ = everywhere in spacetime. In analogy with this fact we define a gauge 
fixing adapted to the weakly isolated horizon by setting: 

£iAf, \a= (77) 



This choice of the gauge is only formally analogous to the condition ( |74[ ) imposed on co (i.e. £iuj = 0) 
because this latter is a restriction on the form of the gravitational field, while the gauge fixing can 
be imposed without constraining the electromagnetic field itself. We can now define the horizon 
electrostatic potential as: 

^(i) = -l^A^ (78) 
and the gauge fixing condition ensures that is constant on the horizon 

£iA^ |a= ^ (i($(i)) |a= (79) 



^ We will not use the above definition in the sequel, but we report it for completeness: 

Definition 4 A weakly isolated horizon becomes an isolated horizon (A, [/]) if it is equipped with an equivalence class 
of null normals to it satisfing the equation: 

[£uD]V = G (76) 

for any vector field V tangential to A and I £ [I], where D denotes the covariant derivative with respect to the metric 
on A. 

We stress however that properties holding for weakly isolated horizons are also true for isolated horizons. In fact it is 
easy to show that every isolated horizon is a weakly isolated horizon (see for details [^). 
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This property, together with the zero principle and the boundary conditions which define the ge- 
ometry of a rigidly rotating horizon, wih ensure that the parameters ^(/), ^(^i), which appear in 
the first principle of thermodynamics for rigidly rotating horizons, are constant on A. 



5.2 Energy for rigidly rotating horizon 



The formalism developed for calculating the variation of energy of the gravitational and the electro- 
magnetic field naturally applies to define the variation of energy of a region of spacetime bounded 
by an isolated horizon. 

We do not mind in the sequel what happens on the outer boundary B; this case has been deeply 
examined in [28| to which we refer the reader for details. We rewrite the expression for energy vari- 
ation in the case of Einstein-Maxwell theory, using a 2-surface with no a priori defined geometry. 
Gluing together formula (|32|) and formula (ra), we find: 



6EiLEM,C) = [ \N6x{V^e)-N''6xiV^fa) + 

JBt ^ 



^-s>^^5xa,.}d'x + 



(80) 



Bt 



[£dV^)5x{9) - 6x{V^)£d9)]d'x 



Bt 



N^6x{E^V^)d? 



Bt 



N''5x{E^KV^)(fx + 



Bt 



N[u^npe'''^'^^'B^,5xAa\^(fx 



This formula expresses the variation of energy for the Einstein-Maxwell theory evaluted on a generic 
surface Bt. From now we shall specialize this formula to the case of rigidly rotating horizons, i.e. 
to the case of Bt = At. We stress that in that case we chose a null evolution of the boundary and 
this implies that the boost velocity, defined in ([27|), is u = 1. 
Using the boost relations (see |l^), we obtain: 

li-flC + jvl] 



— 



where we have defined l^y = —a^UuV a^p, so that / = l^yO^^^ and = rf^V . 
After a long calculation (details can be found in |1C], |28] for the gravitational part, while for the 
electromagnetic part we refer to formula (^)), the equation which defines the variation of energy 
on a generic null surface can be rewritten, starting from (|80| ) , as: 



6Ea{Lem,0 



N5x{^9^i))-N^5x{V^^ 
2 



+ 



SA^''dx(Tttu^(fx + 



At 



N<5^i)6x{E^V^)d^x + - I N''Sx{E^AaV^)d 



(81) 



At 



Nr6xA^V^d^x+ I {{£^pf)6xV^ - {£^V^)SxPr}d- 



where we have set: 



Q/3 



kHf) 
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The acceleration A;^/), in accordance with definition ([7^), is expressed by: 



k^i) = fin'^a^ - u'^ho] + k£i{pf) (82) 



HO 



and the natural connection 1-form tD^ = CTa^fs, in accordance with definition (72), turns out to be: 

OJa = O-a^P = CFa'n'y pU,^ + kdaP^ = -kja + kd^P^ (83) 

To perform the calculations from ( |80| ) to (^ij) we also made use of the following relations, which 
can be easily obtained from (|6^) and (|69[): 



Up to now formula ( |8l|) holds true for any cross section of a generic null hypersurface A. In the 
case of a rigidly rotating horizon it is possible to obtain explicitly a first principle of thermodynamics 
if we recall the more relevant geometric properties we explained in section (5.1), which hold true 
on that particular null surface: 



^(0 |A = 







• k{i^ and <!>(/) are constant on A 

• //. = F^pl^a^^ = and = F^pl'^a^^ = 

• ^" = = - on A 

• (/? is a symmetry for both the 3-metric and the electromagnetic field 

With these assumptions, which follow directly from the definition of rigidly rotating horizon, it is 
straightforward to rewrite (|8^) as: 

5E^{Lem,0 = ^^xAa + <^ii)6xQA + ^(i)SxJa (84) 
where we define the area, the charge and the angular momentum of the horizon as follows: 

Qa = -f J^^Ex^Jad'^x 

= -\ Ia, ^"('^^ + 2E±A^)^d^x 

and owing to the definition of rigidly rotating horizon the quantities and Qa are independent 
on time, i.e. on the slice At we are integrating on. The definition of rigidly rotating horizon ensures 
that also Ja is independent on the slice At we are integrating on; see ||5| and for a detailed 
discussion. 

The integrability of this expression depends on the calibration of the parameters (namely the "tem- 
perature" = the angular velocity r2(;) and the electrostatic potential ^(i))- This argument are 
extensively treated in other papers and we refer for details to Q, [^. We just want to point out that 
the Noether covariant approach together with a proper analysis of boundary terms in the definition 
of energy allows us to reproduce the first law of thermodynamics also in the case of rigidly rotating 
horizons ( p^ ) in exactly the same way as it was formulated in |Q. 
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6 Conclusions 



Our analysis allowed us to formulate a first principle for rigidly rotating horizons analogous to the 
ones given in |^] and but the approach and the formalism used here are completely different. The 
approach used in |^] is metric-affine, based on tetrad gravity, and uses a symplectic type analysis 
to define the conserved quantities. On the other hand, Booth uses in Q a purely metric approach 
based on a canonical decomposition of a "modified" trace-K action functional. Indeed he needs to 
add suitable boundary terms in the action functional to correct the conserved quantities. 
We worked instead in a purely metric framework and we considered a definition of the conserved 
quantities which naturally arises from the Noether theorem. Closed in spirit with the methods em- 
ployed in the covariant ADM approach (|^3[, |^^) and the Regge-Teitelboim analysis of boundary 
terms in the Hamiltonian variation, we are able to correctly define the variation of the Hamiltonian 
and the variation of energy by pushing all the boundary terms arising in the variation of the Noether 
current into the definition itself. 

This allowed us to treat succesfully Einstein-Maxwell theory, where Electromagnetism is considered 
as a natural theory. The framework allows us to define the theory from a geometric point of view, 
without a priori assumptions on the magnetic field and on the existence of a global potential. En- 
ergy definitions arising from the imposed boundary conditions have a nice physical interpretation. 
Finally this theory find applications in defining a first law for rigidly rotating horizons, which have 
recently assumed importance in trying to generalize the quasi static treatments of thermodynamics 
of black holes. These surfaces allow us to define the first principle also for non-static solutions, that 
admit radiation near the horizon; quasilocal formalism is indispensable in this context. 
We have treated here the case of minimal interaction between the gravitational field and the elec- 
tromagnetic radiation using Einstein-Maxwell Lagrangian. It is straightforward to include into the 
theory a dilatonic coupling with the Maxwell Lagrangian and generalize the framework to include 
also more general Yang-Mills fields. This arguments will be treated in a forthcoming paper IQ. 
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